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^ • ABSTRACT 

We present high- resolution (i? ~ 18, 000), high signal-to- noise 2 /im spectra of two luminous, 
I X-ray flaring Class I protostars in the p Ophiuchi cloud acquired with the NIRSPEC spectrograph 

■ of the Keck II telescope. We present the first spectrum of a highly veiled, strongly accreting 
p\J I protostar which shows photospheric absorption features and demonstrates the stellar nature of 

■ its central core. We find the spectrum of the luminous {Lt,oi = 10 Lq) protostellar source, 
(-H ^ YLW 15, to be stellar-like with numerous atomic and molecular absorption features, indicative 
pi i' of a K5 IV/V spectral type and a continuum veiling — 3.0. Its derived stellar luminosity 

(3 Lq) and stellar radius (3.1 Rq) are consistent with those of a 0.5 Mq pre-main-sequence star. 
5^ ' However, 70% of its bolometric luminosity is due to mass accretion, whose rate we estimate to 

^1 be 1.6 X 1Q~^Mq yr^^ onto the protostellar core. We determine that excess infrared emission 

' produced by the circumstellar accretion disk, the inner infalling envelope, and accretion shocks 

at the surface of the stellar core of YLW 15 all contribute significantly to its near-IR continuum 
veiling. Its projected rotation velocity v sin i = 50 km s^^ is comparable to those of flat-spectrum 
. , protostars but considerably higher than those of classical T Tauri stars in the p Oph cloud. The 

■ protostar may be magnetically coupled to its circumstellar disk at a radius of 2 i?, . It is also 
plausible that this protostar can shed over half its angular momentum and evolve into a more 
slowly rotating classical T Tauri star by remaining coupled to its circumstellar disk (at increasing 
radius) as its accretion rate drops by an order of magnitude during the rapid transition between 
the Class I and Class II phases of evolution. 

The spectrum of WL 6 does not show any photospheric absorption features, and we estimate 
that its continuum veiling is rfc > 4.6. Its low bolometric luminosity (2 Lq) and high veiling 
dictate that its central protostar is very low mass, M ^ 0.1 Mq. 

Subject headings: stars:atmospheres, formation, and rotation — infrared: stars — tech- 
niques :spectroscopic 
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1. Introduction 

The discovery of accreting protostars has been 
a triumph for understanding fow-mass star forma- 
tion. Infrared and sub-miUimeter data have re- 
vealed that these objects are deeply embedded, 
cold, and luminous. The spectral energy distri- 
butions (SEDs) of such Class I protostars grad- 
ually rise from the near-IR, peak in the far-IR 
(A ~ 100/um), and fall off in the sub-mm and mm 
wavelength regimes (e.g., Adams, Lada, & Shu 
1987). Observations and models have indicated 
that protostars consist of three primary compo- 
nents: a massive infalling envelope, a circumstellar 
accretion disk, and an embryonic stellar core. 

With a size of ^ 10'' AU, the infalling envelope 
is the largest component of a protostellar object. 
It extends to within a few AU of the embryonic 
stellar core and absorbs and reprocesses all the vis- 
ible (and UV) as well as much of the near-infrared 
radiation emitted by the central protostar. The 
infalling envelope contributes the bulk of the ob- 
served far-infrared to millimeter-wave luminosity 
and a significant portion of the detected near- and 
mid-infrared emission. 

The circumstellar accretion disk can be as much 
as a few hundred AU in size and it contributes a 
significant portion of the observed mid- and near- 
infrared luminosity of the protostar. During the 
protostellar phase of evolution, the embryonic stel- 
lar core at its center grows primarily through ac- 
cretion of material from the circumstellar disk, 
which itself acquires its mass directly from the in- 
falling envelope. 

The embryonic stellar core at the center of the 
protostar is thought to be a stellar-like object a 
few solar radii in extent (e.g., Stabler, Shu, & 
Taam 1980). It is the source of the bulk of the 
protostellar luminosity which is generated by both 
a strong accretion shock at the stellar surface and 
the radiation of internal stellar energy. This inter- 
nal energy is provided by its self-gravity as well as 
the release of gravitational potential energy during 
the protostar's assembly from accreting and in- 
falling material. The core's accretion-dominated 
luminosity is radiated primarily at UV and vi- 
sual wavelengths, but most, if not all of this ra- 
diation is absorbed by dust in the infalling enve- 
lope which renders the central protostar invisible 
at these wavelengths. As a result, little is known 



about the nature of the embryonic stellar cores at 
the hearts of protostellar objects. 

The observational technique of high-resolution 
near-IR spectroscopy has the potential to reveal 
the nature of these heavily embedded central ob- 
jects. Indeed, we have recently been able to 
demonstrate the stellar nature of the central ob- 
jects in flat-spectrum protostars and thus discern 
their spectral types, effective temperatures, rota- 
tion vc;locities, luminosity classes, and masses us- 
ing this technique (Greene & Lada 1997, hereafter 
Paper I). These objects are in a later evolutionary 
state than Class I protostars; they are still sur- 
rounded by substantial circumstellar disks and en- 
velopes but their accretion luminosities (and rates) 
have decreased below Class I levels. However, our 
previous high-resolution IR spectroscopic study of 
Class I young stellar objects (YSOs) did not re- 
veal any spectral absorption features in these pro- 
tostars (Greene & Lada 2000, hereafter Paper II). 
Luhman & Rieke (1999) did detect absorption fea- 
tures in the moderate-resolution near-IR spectrum 
of at least one Class I YSO (YLW 16A / IRS 44), 
but the moderate continuum veiling of that object 
(rk = 1) suggests that it is intrinsically similar to a 
flat-spectrum YSO. Also, Kenyon et al. (1998) de- 
tected TiO absorptions in moderate-resolution op- 
tical spectra of three Class I YSOs in the Tau-Aur 
dark clouds. However, these are all low luminosity 
objects, ~ Lhoi < 0.5 Lq, indicating relatively 
low levels of accretion. Thus there is very little or 
no existing information concerning the nature and 
basic physical characteristics of the central stel- 
lar cores of protostars undergoing significant mass 
accretion. 

Observational determination of the effective 
temperatiircs and radii (or surface gravities) of 
these deeply embedded central objects woiild pro- 
vide critical constraints for theoretical protostellar 
models which predict the values of these parame- 
ters as a protostar evolves (Stabler, Shu, & Taam 
1980). Empirical knowledge concerning the ini- 
tial angular momentum of these objects as well as 
the evolution of angular momentum in the earliest 
phases of young stellar evolution is also highly de- 
sirable to test protostellar theory. Many classical 
T Tauri stars (CTTSs) embedded in the p Oph 
cloud and elsewhere are known to rotate slowly, 
V sin i < 20 km s~^. However, we have found that 
flat-spectrum protostars in this cloud have similar 
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radii but rotate much more quickly, v sin z ^ 40 
km (Paper I). Do heavily-accreting Class I 
protostars in Ophiuchus rotate as fast or faster? 
If so, how is stellar angular momentum reduced so 
quickly by the classical T Tauri pre-main-sequence 
(PMS) evolutionary phase? 

We have undertaken a new, very sensitive near- 
IR spectroscopic study of Class I protostars with 
the 10 m Keck II telescope in order to address 
these issues. The Keck observations are consid- 
erably more sensitive than our previous deep in- 
frared spectral survey made with the 3-m NASA 
IRTF. In this paper we report our findings for two 
of these objects in the p Oph dark cloud, YLW 15 
(IRS 43) and WL 6. These are 2 of the 3 p Oph 
Class I protostars detected in hard X-rays by the 
ASCA satellite, and both showed strong, variable, 
and flaring emission (Kamata et al. 1997). We 
describe our data acquisition and reduction in §2, 
analyze the spectra in §3, discuss the results in §4, 
and summarize our conclusions in §5. 

2. Observations and Data Reduction 

Near-IR spectra of the p Oph protostar YLW 
15 were acquired on 2001 July 7 UT, and the pro- 
tostar WL 6 was observed on 2000 May 30 UT. 
MK spectral standards were also observed 2000 
May 29 - 30 UT. All data were acquired with the 
10-m Keck II telescope on Mauna Kea, Hawaii, 
using the NIRSPEC facility multi-order cryogenic 
echelle spectrograph (McLean et al. 1998). Spec- 
tra were acquired with a 0'.'58 (4 pixel) wide slit, 
providing spectroscopic resolution R = X/5X = 
18,000 (16.7 km s^^). The plate scale was 0'.'20 
pixel"^ along the 12" slit length, and the see- 
ing was 0'.'5 - O'.'e. The NIRSPEC gratings were 
oriented to allow both the 2.21 pm Na and the 
2.30 pm CO band head regions to fall onto the 
instrument's 1024 x 1024 pixel InSb detector ar- 
ray, and its NIRSPEC-7 blocking filter was used. 
NIRSPEC was configured to acquire simultane- 
ously multiple cross-dispersed echelle orders 31 
- 36 (2.08 2.45 pm, non-continuous) in 2000 
May and orders 32 - 38 (1.97 - 2.38 ^m, non- 
continuous) on 2001 July 7. Each order had an 
observed spectral range AA ~ A/67 {Av ~ 4450 
km s~^). 

The internal instrument slit rotator was used to 
maintain a fixed position angle on the sky when 



observing WL 6 and the MK standards on 2000 
May 29 - 30. The slit was held physically sta- 
tionary and thus allowed to rotate on the sky (as 
the non-equatorially-mounted telescope tracked) 
when observing on 2001 July 7 (YLW 15). Data 
were acquired in pairs of exposures of durations 
from less than 1 second (for giant MK standards) 
to up to 600 s (for the YSOs) each, with the tele- 
scope nodded 6" along the slit between frames 
so that object spectra were acquired in all expo- 
sures. WL 6 was observed for a total of 30.0 min- 
utes, while YLW 15 was observed for 73.3 minutes. 
Early-type (B7 - A2) dwarfs were observed for tel- 
luric correction of the MK standard stellar spec- 
tra. HR 6054 (B7 IV) and HR 6070 (AO V) were 
observed for telluric correction of the WL 6 and 
YLW 15 spectra, respectively. The telescope was 
automatically guided with frequent images from 
the NIRSPEC internal SCAM IR camera during 
all exposures of more than several seconds dura- 
tion. Spectra of the internal NIRSPEC continuum 
lamp were taken for flat fields, and exposures of 
the its Ar, Ne, Kr, and Xe lamps were used for 
wavelength calibrations. 

All data were reduced with IRAF. First, object 
and sky frames were differenced and then divided 
by flat fields. Next, bad pixels were fixed via in- 
terpolation, and spectra were extracted with the 
APALL task. Extracted spectra were typically 4 
pixels (0'.'6) wide along the slit (spatial) direction 
at their half-intensity points. Spectra were wave- 
length calibrated using low-order fits to lines in the 
arc lamp exposures, and spectra at each slit po- 
sition of each object were co- added. Instrumental 
and atmospheric features were removed by divid- 
ing wavelength-calibrated object spectra by spec- 
tra of early-type stars observed at similar airmass 
at each slit position. Combined spectra were pro- 
duced by summing the spectra of both slit posi- 
tions for each object and then shifting the spectra 
in wavelength to correct for radial velocities. The 
standard spectra were shifted so that their lines 
matched laboratory wavelengths, while the YSO 
spectra were corrected for all Earth and solar sys- 
tem radial motions relative to the local standard of 
rest. The YSO spectra were not further corrected 
for the radial velocity of the p Oph cloud. Fi- 
nally, YSO spectra were multiplied by an artificial 
spectrum of a T=10,000K black body in order to 
restore their true continuum slopes (although this 
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was a small correction of under 10% per order). 

3. Data Analysis and Results 
3.1. Spectra 

The spectra of several MK standard stars are 
shown in Figure 1. The three NIRSPEC orders 
with the most numerous absorption features are 
displayed. These features are similar to ones seen 
by Wallace & Hinkle (1996) in stars of similar 
spectral type, although we present a finer grid 
of spectra in the K2V M2V range. Spectra of 
these stars arc dominated by strong features of 
Si (2.0923 Aim), Mg (2.1066 ^im), and Al (2.1099 
lim) in the shortest wavelength order, with ratios 
Al/Mg and Al/Si being very sensitive to spectral 
type. Some of the strongest lines in the next or- 
der are Na (2.2062, 2.2090 ^^m), Sc (2.2058, 2.2071 
Hm), and Ti (primarily 2.2217 and 2.2240 jim). 
The Ti lines are very good indicators of spectral 
type in this range, growing stronger from earlier 
to later types. The 2.2069 /im Si hue is also a 
good indicator of spectral type, weakening from 
early to later types. Therefore the ratio of Ti/Si 
is a particularly good spectral type diagnostic in 
this order. The strongest features in the longest 
wavelength order are Mg (2.2814 /im) and the CO 
V — — 2 band head and rotation-vibration lines 
(2.2935 nm and long- wards). The Mg line becomes 
significantly weaker in later spectral types. All 
transitions identified above and hereafter (and in 
Fig. 1) occur in neutral species. 

The two giant spectra (K3 III and M2 III) 
shown in Figure 1 have many of the same lines 
as the dwarfs, but there are significant differences. 
The Si and Mg lines in the shortest order do 
not appear to change monotonically with spectral 
type, but the lines in the next order are much 
more telling. The Sc (2.2058, 2.2071 /xm) lines 
are very strong (absolutely and relative to Na) in 
the M2III giant, while the Ti lines are very strong 
in both giants. The Na lines (2.2062 and 2.2090 
/xm) are less broad in the giants because they are 
formed under conditions of lower pressure than 
in dwarfs. The Mg (2.2814 iim) line is also less 
broad in the giants and does not vary monotoni- 
cally with spectral type. The 2.2935 /im CO band 
head and nearby rotation-vibration lines are very 
strong in the giants. Numerous other low-strength 
absorptions are seen in all three orders of the gi- 
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Fig. 1.— Spectra of MK standards. Three NIR- 
SPEC orders are shown, and strong absorption 
features are identified in each order. Each of the 
three orders contain a number of lines whose ra- 
tios are useful for evaluating the spectral types 
of protostars independent of their veiling. Each 
spectrum has been normalized to a value of 1.0 
and they are offset vertically by a constant in each 
order for separation. 



4 



ant star spectra, indicating the presence of various 
molecules (see also Wallace & Hinkle 1996). 

Three NIRSPEC spectral orders of the YSOs 
WL 6 and YLW 15 are shown in Figure 2. The WL 
6 data have signal-to-noise of 100 - 170 at most 
wavelengths, and the YLW 15 data have signal- 
to-noise 240 - 300. Both objects have strong con- 
tinua rising towards longer wavelengths. YLW 15 
has numerous absorption features which are weak 
and moderately broad, while any absorption fea- 
tures in the spectra of WL 6 are either weaker or 
non-existent. Both sources exhibit strong v = 1 
- S(0) H2 emission at 2.2233 ^m. Interestingly, 
the strength of this emission appears to be corre- 
lated with infrared excess and veiling in these ob- 
jects. WL 6 has weaker (absent) absorption lines, 
a stronger red continuum, and stronger H2 emis- 
sion than YLW 15. 

3.2. Veiling and Rotation Analysis 

The spectra of YLW 15 and WL 6 are com- 
pared to those of the K4 dwarf HR 5568 in Fig- 
ure 3. In this figure, the spectra of the YSOs are 
processed beyond that shown in Figure 2 in order 
to facilitate this comparison. First, their contin- 
uum slopes were removed by normalization. Next, 
the relative strengths of their absorption lines were 
magnified by subtracting off a constant continuum 
value of 0.75, 75% of the normalized continuum. 
This specific constant was chosen so that the ab- 
sorption lines would have nearly the same equiv- 
alent widths in both YLW 15 and HR 5568. The 
spectrum of HR 5568 was modified by convolv- 
ing it with a limb-darkened stellar rotation profile 
with velocity v sin i ~ 50 km s~^. This was re- 
quired to match the velocity widths of its absorp- 
tion features to those of YLW 15. 

Figure 3 shows that the spectrum of YLW 15 
is very similar to that of HR 5568. Nearly all 
absorption features have nearly identical relative 
strengths. Comparing relative line strengths is 
most useful because ratios (but not absolute equiv- 
alent widths) are insensitive to continuum veiling. 
For example, the ratios of their Mg (2.1066 fim) 
and Al (2.1099 /xm) equivalent widths are very 
similar, as are the ratios of their lines in the Na and 
Sc region (2.2053 2.2096 /im). The Hues in the 
Na and Sc region also have the same strengths rel- 
ative to the Fe and Sc lines in the 2.2263 - 2.2273 
lim region. However, the spectra of the two ob- 
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Fig. 2. — Spectra of Class I protostars YLW 15 
and WL 6. The absorption lines are either weak 
(YLW 15) or undetected (WL 6) due to the large 
continuum veilings (r^ > 3) of these objects. Each 
spectrum has been normalized to a value of 1.0 
(over all 3 orders), and WL 6 has been offset ver- 
tically by 0.2 for in each order for separation. 
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Fig. 3. — Spectra of YLW 15 and WL 6 compared 
to HR 5568 (K4 V). The spectra of the protostars 
have been normalized to remove their continuum 
slopes, and 75% of their continuum levels have 
been subtracted to remove the continuum veiling 
amount rj. = 3. Each spectrum has been nor- 
malized to a value of 1.0. The spectra are offset 
vertically by a constant in each order for separa- 
tion. 



jects do not match exactly. The 2.2280 and 2.2317 
Hm lines of Ti are stronger relative to the 2.2053 

- 2.2096 /Ltm Na, Sc, and Si lines in YLW 15 than 
in HR 5568. Likewise, the 2.2814 ^m Mg line is 
weaker relative to the CO band head (or the 2.2053 

- 2.2096 Atm Na, Sc, and Si lines) in YLW 15 than 
in HR 5568. Both of these differences suggest that 
YLW 15 is of slightly later spectral type than HR 
5568. It is clear from Figure 1 that YLW 15 is not 
as late as Gl 338B (K7 V) and that its spectrum 
is much more dwarf-like than giant-like (see also 
the discussion of giant spectra in §3.1). We assign 
the spectral type K5 IV/V to YLW 15. Differ- 
entiating between luminosity class IV and V will 
probably require comparisons to synthetic model 
spectra given the rarity of well-classified luminos- 
ity class IV stars. 

The continuum veiling of YLW 15 can be mea- 
sured now that its spectral type has been de- 
termined. Continuum veiling at wavelength A 
is defined as the value of continuum excess ex- 
pressed as a fraction of the stellar photospheric 
flux, rx = Fx^ex/Fx*- We have estimated the K- 
band continuum veiling of YLW 15, rfc, by com- 
paring the equivalent width of its K-band spectral 
features to those of MK standards. The above con- 
tinuum veiling equation can be applied to K-band 
wavelengths and written in terms of the intrinsic 
{EW) and observed (EW) equivalent widths of 
spectral features there, EW' = EW/ (1 -I- r^). 

We have evaluated of YLW 15 by compar- 
ing the measured equivalent widths of its 2.2053 

- 2.2096 ^m Na, Sc, and Si lines (from its spec- 
trum in Figure 2) to those of the K4 V star HR 
5568 and the K7 V star Gl 338B. The measured 
equivalent widths of these features are 0.8, 2.8, and 
3.5 A for similarly-reduced spectra of YLW 15, 
HR 5568, and Gl 338B, respectively. The intrinsic 
equivalent widths of these features is expected to 
be approximately 3.2 A for YLW 15, intermedi- 
ate between HR 5568 and Gl 338B. Substituting 
these numbers into the above equation yields rk 
= 3.0. We estimate that this value is accurate to 
about 25% given the uncertainties in the measured 
equivalent widths and the spectral type determi- 
nation of YLW 15. This veiling implies that 75% 
of the iC-band flux of YLW 15 arises from (cir- 
cumstellar) continuum excess, and 25% of its fltix 
arises from its K5 stellar photosphere. Therefore 
the stellar spectrum of this object can be approx- 
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imated by subtracting 75% of the continuum level 
at if -band wavelengths. This - and normaliza- 
tion by its continuum - was done to produce the 
spectrum of YLW 15 shown in Fig. 3. 

The derived veiling rk = 3.0 of YLW 15 is con- 
siderably less than the rfc > 6 which we estimated 
for this source in Paper II. There may be several 
reasons for this discrepancy. First, we have re- 
examined the older, lower quality data of YLW 15 
and have concluded that CO band head absorp- 
tion may be present at low confidence in its spec- 
trum (see Figure 2 of Paper II). The presence of 
such absorption would have lowered our estimate 
of the amount of veiling for this source. Second, 
the veiling in the protostar may have decreased 
during the 4 years which elapsed between our two 
observations. Luhman & Rieke (1999) found that 
p Oph sources can undergo large variations in veil- 
ing, and they also derive > 1 from their 2 /xm 
spectrum of YLW 15. Finally, in Paper II we as- 
sumed that YLW 15 was spectral type MO when 
we estimated its veiling. Its actual spectral type is 
K5 (see above), which is expected to have intrin- 
sically weaker CO absorptions (analyzed in paper 
II) than would a MO star. This would have led to 
an overestimate of its veiling. 

We have determined that the projected rota- 
tion velocity of the YLW 15 central protostar is v 
sin i = 50 ± 5 km a~^. This was measured by fit- 
ting its CO band head spectrum (Fig. 3c) to the 
spectrum of HR 5568 (K4 V) which had been con- 
volved with stellar rotation profiles over a range of 
velocities. In order to confirm that this value was 
not influenced by any slight rotation of HR 5568 
(whose rotation is unresolved in our spectra), we 
also fit the spectrum of YLW 15 with that of a 
slowly-rotating K giant star which had been con- 
volved with stellar rotation profiles over a range 
of velocities. This also resulted in a projected ro- 
tation velocity of 50 km s~^. 

The spectra of WL 6 (Figs. 2 and 3) do not 
show any obvious absorption features, so only a 
lower limit to its veiling can be calculated. If the 
central protostar of WL 6 has the same spectral 
type and v sin i rotation velocity as YLW 15, then 
the signal-to-noise of its spectrum (^ 120 in the 
2.2053 - 2.2096 /zm region) dictates that WL 6 
must have > 4.6. This is consistent with the 
lower limit determined in Paper II, rk > 3.8 — 5.5. 
This veiling analysis was based on maximum line 



depths (and not equivalent widths), so we expect 
the veiling of this object to be changed little if it is 
of later spectral type than YLW 15. Our analysis 
of dwarf standards indicates that although equiv- 
alent widths of the 2.2053 - 2.2096 /xm features in- 
crease in later types (particularly in the M2 - M6 
range), this is mostly caused by increasing intrin- 
sic line widths (due to pressure broadening), with 
the peak line depths relatively unchanged. How- 
ever, less veiling is required if WL 6 is rotating 
more rapidly than YLW 15, and more is required 
if it is rotating more slowly. Rotation can affect 
the measured line depths and veilings by about a 
factor of 2 over the range v sin i = - 180 km s~^ 
(see Paper II). 

4. Discussion 

We have determined that the Class I YSOs WL 
6 and YLW 15 have strong continuum veilings, 
and both have near-IR H2 emission. Absorption 
lines in the spectrum of YLW 15 have allowed us to 
measure quantitatively its spectral type, rk contin- 
uum veiling, and v sin i rotation velocity. We now 
examine and discuss the physical implications of 
these results given other information known about 
these YSOs. 

4.1. YLW 15 

The SED of YLW 15 rises steeply from the 
near- to the far-IR (2 60 /im) with a bolomet- 
ric luminosity of approximately 10 L© (Wilking 
et al. 1989). Its 2-14 ^m spectral index (slope 
of its SED) is 1.0 (Bontemps et al. 2001). Andre 
& Montmerle (1994) found that this source also 
has significant A = 1.3 mm spatially-extended flux 
(100 mjy), which they interpreted as arising from 
its circumstellar envelope, confirming its Class I 
nature. Bontemps et al. (1996) found that YLW 
15 has a compact molecular outflow which is seen 
nearly pole-on. The coexistence of molecular out- 
flow and near-IR H2 emission (Figs. 2 and 3) in 
this object is consistent with the correlation of 
these phenomena found by Greene & Lada (1996). 
GY 263 is a fainter near-IR source {6K = 2.5 mag; 
Barsony et al. 1997) 6" (~ 900 AU projected dis- 
tance) NW of YLW 15. 

However, Girart, Rodriguez, & Curiel (2000) 
recently found that YLW 15 itself is binary in their 
A = 3.6 cm VLA map. 2MASS (second incremen- 
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tal release) Point Source Catalog coordinates show 
that the 2 /zm source is essentially coincident with 
VLA 2, the weaker, point-like radio source. VLA 
1 is located approximately 0'.'6 100 AU) to the 
NW and is spatially extended (~ 1") SW - NE. 
Girart, Rodriguez, & Curiel (2000) interpret the 
VLA 1 emission to be thermal free-free, arising in 
a protostellar jet. HST / NICMOS near-IR im- 
ages (FllOW and F160W bands) show that the 
central soiirce is single but embedded in extended 
nebulosity (Allen et al. 2002). Hereafter we as- 
sume that there is a central protostellar engine for 
YLW 15, and it is located at the VLA 2 / near-IR 
position. The VLA 1 source is likely to be emis- 
sion from a jet produced by the central protostar 
and not a (proto-) stellar source itself. 

One of the most intriguing aspects of YLW 15 
is its X-ray emission. This was one of the very few 
Class I protostars which was detected by ROSAT 
(Casanova et al. 1995). Subsequent ROSAT HRI 
(Grosso et al. 1997), ASCA (Tsuboi et al. 2000), 
and Chandra (Imanishi, Koyama, & Tsuboi 2001) 
observations have all shown that its X-ray emis- 
sion has frequent strong flares (~ 10^^ erg s^^; see 
Grosso et al. 1997; Imanishi, Koyama, & Tsuboi 
2001). The ASCA and Chandra observations have 
shown that its quiescent and flaring emissions con- 
sist of mostly hard X-rays. The ASCA observation 
captured 3 quasi-periodic flares with a ~ 20 h pe- 
riod which spanned the exposure (Tsuboi et al. 
2000), but the Chandra exposure captured only a 
single flare ('^ 3 hr duration) in a 27.8 h exposure. 
Montmerle et al. (2000) have interpreted this flar- 
ing as arising from the interaction of the magnetic 
fleld of the rapidly rotating protostar (P ~ 20 h, 
near break up) and its slower rotating circumstel- 
lar disk. 

Our near-IR spectrum of YLW 15 provides fur- 
ther information on the nature of this object. We 
have assigned it a spectral type of K5 IV/V (see 
§3.2), implying an effective temperature Teff — 
4300 K (Tokunaga 2000). Its near-IR spectrum 
(Fig. 3) shows no evidence for binarity; it shows 
no features which could arise in a companion of 
different spectral type. We now derive additional 
stellar and system parameters for this object. 

4.1.1. Stellar Parameters 

We used new near-simultaneous near-IR pho- 
tometry (acquired 2001 July 8; see Haisch et al. 



2002) to estimate the extinction towards the cen- 
tral protostar of YLW 15. The observed CIT 
magnitudes of this object were H = 12.53 and 
K = 9.35 (undetected at J), and we de- reddened 
the H — K color to determine that the extinc- 
tion to the central protostar is ~ 40 mag. 
This value is intermediate between the lower limit 
[Ay = 34 mag), resulting from de-reddening to 
the maximum CTTS locus color [H - K)o = 1.0 
mag (sec Meyer, Calvct, & Hillenbrand 1997), and 
the upper limit {Ay = 47 mag), resulting from 
de-reddening to the intrinsic stellar photospheric 
color {H - K)o = 0.11. YLW 15 cannot be ac- 
curately de-reddened to this stellar photospheric 
color because the de-reddened colors of Class I 
YSOs do not intersect stellar loci in near-IR color- 
color diagrams. This is due to the presence of both 
scattered light and infrared excess emission from 
the close-in circumstellar material contained in the 
surrounding accretion disk and infalling envelope 
(see Lada & Adams 1992). Our adopted extinc- 
tion {Ay ~ 40 mag) is consistent with the Ay > 31 
mag value derived by Grosso et al. (1997) using 
a similar technique. Our extinction value, along 
with the measured near-IR excess = 3.0, the 
IR extinction law, dwarf bolometric corrections, 
and a distance modulus of 6.0 yields an estimated 
stellar luminosity of 3.0 Lq. The 4300 K effec- 
tive temperature of the protostar dictates that its 
stellar radius is 3.1 Rq. Although there are con- 
siderable uncertainties in extinction (and scatter- 
ing), the derived protostellar properties are not 
greatly affected. Varying the Ay between its lower 
and upper limits changes the derived stellar lumi- 
nosity by less than a factor of 2 and the derived 
stellar radius by less than 50%. The derived lu- 
minosity and radius are very similar to those pre- 
dicted for a young PMS star of this effective tem- 
perature which is still accreting mass at a high 
rate (M ~ IO^^Mq yr^^; see Siess, Forestini, & 
Bertout 1999; Tout, Livio, & Bonneh 1999). These 
and other modern PMS models indicate that the 
mass of YLW 15 is approximately 0.5 M©, as de- 
termined from its effective temperature and stel- 
lar luminosity (D'Antona & Mazzitelli 1997; Siess, 
Forestini, & Bertout 1999; Siess, Dufour, & Fores- 
tini 2000). Taken with its derived radius, this dic- 
tates that the surface gravity of YLW 15 is Log 
5 ~ 3.1 (cgs units), within the range expected for 
a very young K5 PMS star. 
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4.1.2. Accretion 

The derived stellar luminosity is about 30% of 
YLW 15's bolometric luminosity (Wilking et al. 
1989), implying that 70% of the total luminosity 
of this protostar is due to mass accretion and not 
radiation of the protostar's internal energy. This 
accretion luminosity is likely to be composed of a 
component which arises in a circumstcllar accre- 
tion disk as well as a component which is gener- 
ated when accreting matter falls directly onto the 
central protostar. If all matter accretes through 
a circumstellar disk with inner radius Rin, then 
the accretion luminosity arising in the disk as 
GMMi)/2Rin where Md is the mass accretion 
rate through the disk. The accretion luminosity 
produced by matter falling from the disk edge Rin 
onto the star is GMM^/R^ - GMM^JRin where 
M* is the mass accretion rate onto the star. We 
now assume that matter is transported from near 
the disk edge to the stellar surface in magneto- 
spheric accretion columns (Koenigl 1991; Shu et 
al. 1994). In this scenario, matter accretes through 
the disk to a disk coupling radius R^ where it is 
then accreted onto the star along magnetic field 
lines (Shu et al. 1994). This point R^ is slightly 
exterior to Rin, and the total accretion luminosity 
of the system is 

_ gmMd gmja^ _gmJu 

The typical inner hole radius of a CTTS cir- 
cumstellar disk is Rin — 5i?* (Shu et al. 1994; 
Meyer, Calvet, & Hillenbrand 1997). Magneto- 
spheric accretion models predict that the inner 
disk and coupling radii are closer to the central 
star when matter is accreting at higher rates, such 
as those typical for Class I protostars. DQ Her 
stars are also sometimes interpreted as exhibiting 
this behavior (e.g., see Patterson 1994). We as- 
sume M* = 0.7Md from Shu et al. (1994) and 
adopt their equation (2.6b) to express the mag- 
netic coupling radius as 

Rx f B^R^ \ 

where ax is a constant on the order of 1 (Shu 
et al. 1994; Ostriker & Shu 1995). 



We solved equations (1) and (2) simultaneously 
for mass accretion rates and Rx using the previ- 
ously derived values of i?» (3.1 Rq), M* (0.5 Mq), 
ax = 0.92 (Ostriker & Shu 1995), the above rela- 
tion of mass accretion rates, and the assumption 
= 1000 G. This results in a magnetic cou- 
pling radius Rx = 2.1R^ and mass accretion rates 
Md =2.3xlO-6MQyr-iandM* = 1.6x10-^Mq 
yr~^. This magnetic coupling radius is consider- 
ably uncertain (perhaps by 50%) because we have 
assumed a magnetic field strength and have not 
measured a value. The 2.2233 /xm TI line is sen- 
sitive to magnetic (Zecman) broadening, and in- 
deed it has a larger equivalent width (~ 0.4 A) 
in the stellar spectrum of YLW 15 than in either 
the K4 V standard HR 5568 (0.2 A) or the K7 
V standard Gl 338B (0.3 A; see Figs, lb and 3b). 
This suggests the presence of a significant (kG; e.g. 
sec Johns-Krull, Valenti, & Koresko 1999) mag- 
netic field in YLW 15, but the magnitude of this 
field cannot be measured quantitatively from our 
spectrum due to its limited signal-to-noise and the 
presence of the flanking 2.2233 /xm H2 line (Fig. 
3b). 

The derived accretion rates are consistent with 
the accretion rates computed for other Class I 
YSOs in this cloud and are within a factor of 2 
of the spherical accretion rate derived from the 
cloud's gas temperature (M = a?/G; see Adams, 
Lada, & Shu 1987). The adopted 0.5 Mq stellar 
mass of YLW 15 is also consistent with accretion 
at this rate (i.e. M*) over the mean Class I life- 
time in the cloud (Greene et al. 1994; Luhman & 
Rieke 1999). The derived accretion rates and in- 
ner disk radius result in accretion luminosities of 
2.75 Lq arising in the disk and 4.25 Lq arising at 
the stellar surface. Thus the stellar surface accre- 
tion luminosity is approximately 40% greater than 
the luminosity emanating from the Teff = 4300 
K protostellar photosphere. Calvet & GuUbring 
(1998) have shown that the region where an accre- 
tion column meets the surface of a T Tauri star has 
an effective temperature of approximately 8000 K. 
From this we calculate that the derived 4.25 Lq 
stellar surface accretion luminosity arises from an 
area only 12% of the stellar surface. 

4.1.3. Veiling 

Here we investigate the possible sources of ex- 
cess infrared continuum emission and their contri- 
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butions to the overall veiling in the infrared spec- 
trum of YLW 15. We consider three sources for 
the continuum veiling: the accretion luminosity 
generated at the surface of the protostar. the lu- 
minous accretion disk, and the inner infalling en- 
velope. We calculate that the continuum from the 
hot accretion impact sites at the surface of the pro- 
tostar will produce an effective veiling rk = 0.33. 
Although a relatively small fraction of the total 
veiling, this is significantly greater than that pre- 
dicted {rk < 0.05) for accretion impact sites on 
typical class II T Tauri stars whose accretion rates 
arc substantially lower (Calvet & GuUbring 1998). 
Accretion models predict that the disk of YLW 15 
{Ldisk = 2.75^0) should contribute an IR veiling 
Tk ^ 1.3 (see Fig. 9 of Greene & Lada 1996). 
The remainder of the object's observed veiling 
(rfe ~ 1.4 of rfc =3 total) is likely to arise from 
grains in the circumstellar envelope of the proto- 
star which are warmed by the photons produced 
by the stellar photosphere and at the hot sites of 
mass accretion. This implies that the envelope 
contains a significant amount of dust within ~ 0.3 
AU of the central star (see also Greene & Lada 
1996; Calvet, Hartmann, & Strom 1997). 

We note that the envelope veiling is low enough 
so that any absorption lines originating in the cir- 
cumstellar disk should have been seen in the spec- 
trum of YLW 15. However, the (stellar) dwarf- 
like natures of the YLW 15 absorption features 
imply that they are not formed in a circumstellar 
disk, which would show evidence of much lower 
surface gravity (see Paper I for a discussion and 
examples of FU Ori disk features). We calculate 
that the dust grain destruction distance (where 
T < 1800 K) for YLW 15 is approximately 10 
i?* from considering both its protostcUar photo- 
spheric and stellar accretion luminosities. This 
implies that the disk is completely gaseous and 
dust-free between Rin = 2i?* and 10 i?*. The fact 
that we do not observe absorption lines arising in 
this disk region strongly implies that the gas in 
this rc;gion docs not have a siifficdent vertical tem- 
perature gradient to form strong spectral lines. 

4.. 1.4- Rotation 

The V sin i = 50 km s~^ projected rotation 
velocity of YLW 15 is considerably greater than 
typical for many classical T Tauri stars, about 
equal to that found for flat-spectrum protostars 



in the p Oph cloud (see Paper I), and consid- 
erably less than the equatorial rotation velocity 
deduced for this object from its X-ray emission. 
If the 20 h X-ray flare period observed by Mont- 
merle et al. (2000) was caused by rotation as they 
postulated, then assuming a 3.1 Rq protostellar 
radius would imply an equatorial rotation veloc- 
ity of approximately 180 km s~^, approximately 
break up. If YLW 15 is rotating this fast, then 
its rotation axis must be only 16° from our line 
of sight. This is consistent with its molecular out- 
flow being oriented nearly pole-on (Bontemps et 
al. 1996) if its outflow axis is aligned with its ro- 
tation axis. However, subsequent Chandra X-ray 
observations of YLW 15 do not support the inter- 
pretation that it is rotating with a 20 h period. 
These observations did not confirm fiaring with 
a 20 h period; only one flare was detected in a 
10^ s (27.8 h) observation (Imanishi, Koyama, & 
Tsuboi 2001). This weakens the possibility that 
the ASCA X-Ray flare period is physically related 
to the protostellar rotation rate. The inclination 
- and hence the equatorial rotation velocity - of 
YLW 15 is not well constrained. In order to make 
quantitative comparisons, we assume that the ro- 
tation axis of YLW 15 is inclined 57° from our line 
of sight, the most likely inclination angle for a ran- 
dom distribution. This implies that its equatorial 
rotation velocity is approximately 60 km s~^, but 
this is highly uncertain: our newly-determined pa- 
rameters and previous data do not constrain this 
velocity well within the range Vrot = 50 - 180 km 

We estimate that the Keplerian co-rotation ra- 
dius of the YLW 15 circumstellar disk is 2.0 i?* (for 
Vrot = 60 km s~^), nearly equal to the R = 2.1i?* 
magnetic coupling radius. The co-rotation radius 
varies over the range of 2.3 - 1.0 i?* for the range 
of equatorial rotation velocities Vrot = 50 - 170 
km s~^. Due to its large uncertainty (see §4.1.2), 
we cannot determine whether the magnetic cou- 
pling radius is significantly different from the co- 
rotation radius over a large range of inclination 
angles. However, it appears that inclination an- 
gles i ~ 65° (similar to our 57° assumption) are 
most compatible with the Shu et al. (1994) require- 
ment that the co-rotation and magnetic coupling 
radii be essentially equal (provided that we have 
adopted the correct value of i?*). 

If the central protostar and the disk are suf- 
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ficiently magnetically coupled in this region, then 
the disk may be regulating the angular momentum 
(and rotation velocity) of the star (e.g., Koenigl 
1991; Edwards ct al. 1993; Shu ct al. 1994). Al- 
though the observational evidence of such coupling 
in CTTSs has been recently challenged (e.g., see 
Stassun et al. 1999, 2001), the similarity of the 
co-rotation and coupling radii of YLW 15 are con- 
sistent with stellar spin-down and magnetic star- 
disk coupling. Indeed, we calculate that the mag- 
netic braking time scale, r > \ Jt,/{dJ/dt)\ (F. Shu 
priv. communication; eq. 6 of Hartmann 2002), is 
only > 14,000 yr for YLW 15. This is significantly 
shorter than the ~ 10^ yr age of Class I protostars, 
so YLW 15 is likely to have had more than enough 
time to couple magnetically to its disk. 

Class II YSOs and CTTSs in the p Oph cloud 
(and CTTSs in Taurus; see Edwards et al. 1993) 
are observed to rotate considerably more slowly 
{v sin i ~ 12 km s~^; see Paper I), so if YLW 
15 is to evolve into a typical Ophiuchi CTTS it 
will likely lose about 75% of its angular momen- 
tum. This can result via disk braking if the pro- 
tostar remains coupled or recouples to its disk as 
its mass accretion rate drops. If the mass accre- 
tion rate of YLW 15 suddenly dropped to 10"'' 
yr"^, then it would re-couple to its disk in only 
about 150,000 yr (i.e. eq. 6 of Hartmann 2002). 
This coupling would occur at a radius R ~ 4.5i?* 
with co-rotation speed t; sin z ~ 15 km (for 
i = 57°), very similar to the mean rotation velocity 
observed for p Oph Class II YSOs. We note that 
this 150,000 year coupling time is much shorter 
than the mean age of p Oph Class II YSOs (5 - 10 
xlO^ yr), so this can explain why Class II YSOs 
are observed to rotate more slowly than Class I 
and flat-spcctrmn ones. However, more sophisti- 
cated modeling is required to ascertain the exact 
nature of this spin-down behavior. 

4.2. WL 6 

Bontemps et al. (2001) found that WL 6 has 
a spectral index of 0.6 (over A = 2 - 14 /itm), 
less than that of YLW 15. Wilking ct al. (1989) 
also found that the WL 6 spectral index is less 
than YLW 15, but their absolute values and the 
spanned wavelength range are different. WL 6 
has a bolometric luminosity of approximately 2 
L0 (Wilking et al. 1989; Bontemps et al. 2001). 
Neither Andre & Montmerle (1994) or Motte, An- 



dre, & Neri (1998) detected WL 6 at A 1.3 mm 
(peak < 20 mJy in 11" diameter), but both in- 
vestigations regarded it as an embedded protostar 
due to its large spectral index. This source is not 
known to drive a molecular outflow (Bontemps et 
al. 1996), but it does have near-IR H2 emission 
(Figs. 2 and 3). WL 6 source has no known com- 
panions, and it also appeared single at A = 2 /xm 
in the NIRSPEC camera images. 

Kamata et al. (1997) found that WL 6 had 
a hard X-ray spectrum and its flux varied sinu- 
soidally with a 23 h period, which they attributed 
to protostellar rotation. Montmerle et al. (2000) 
also modeled this variation as being due to rota- 
tion, but with a slower period (3.3 d). Imanishi, 
Koyama, & Tsuboi (2001) did not detect this si- 
nusoidal variation, but they did detect two weak 
flares during their 10^ s Chandra exposure. Taken 
with the results of YLW 15, we interpret these dif- 
ferences to indicate that the X-ray variability of 
WL 6 is not directly related to its protostellar ro- 
tation period. However, this cannot be confirmed 
until the rotation signature of a photospheric fea- 
ture (such as 2.29 /xm CO absorption) is observed 
in this object. 

We have used the bolometric luminosity of WL 
6 and the accreting PMS models of Siess, Forestini, 
& Bertout (1999) to bound the stellar parameters 
of this object. The stellar luminosity of WL 6 must 
be less than or equal to its bolometric luminosity, 
2 Lq. The Siess, Forestini, & Bertout (1999) mod- 
els predict that such a PMS star at a very young 
age (~ 10^ yr) must have mass < O.3A'/0 and 
Teff < 4000 K (their Case B). The high veiling of 
WL 6 (see §3.2) implies that its stellar luminosity 
is considerably lower than 2 Lq. We now assume 
that WL 6 has approximately the same ratio of 
stellar to bolometric luminosity of YLW 15 since 
the rfc of WL 6 is somewhat larger and its Teff is 
somewhat lower than YLW 15, with these effects 
tending to cancel each other. This implies that 
the stellar luminosity of WL 6 is ~ 0.6^0, with 
corresponding ~ O.IM0 and Teff < 3160 K. 
This temperature is consistent with that of a M3 - 
M4 PMS star. Our previous models show that it is 
possible to produce significant veiling from only a 
circumstellar accretion disk for a PMS star of this 
mass (see Fig. 9 of Greene & Lada 1996). How- 
ever, it is also plausible and possible that matter 
may also accrete directly from the disk to the star. 
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and some veiling could also be produced by its 
circumstellar envelope (all as in YLW 15). These 
limits on the mass, effective temperature, and lu- 
minosity of WL 6 imply that its mass accretion 
rate is on the order of 10~^ Mq yr~^, from the re- 
lation Lace = Lboi — = GM^M / Rt,. This rate 
is consistent with the typical ~ 10^ yr Class I life- 
time and the ~ 0.1 mass of WL 6. 

5. Summary and Conclusions 

We have presented the first spectrum of a highly 
veiled, strongly accreting Class I protostar (YLW 
15 / IRS 43) which shows late-type photospheric 
absorption features. This spectrum has allowed us 
to diagnose a plethora of details about the stellar, 
accretion, and rotation properties of this proto- 
star and its circumstellar environment. The spec- 
trum of another Class I protostar, WL 6, did not 
show any photospheric absorption features. Both 
of these objects are known to have strong and vari- 
able hard X-ray emission. Our analysis of these 
spectra results in the following conclusions: 

1. We have detected numerous late- type ab- 
sorption features in the photosphere of YLW 15. 
These features include Si (2.0923 jim. and 2.2069 
urn), Mg (2.1066 /xm), Al (2.1099 jim), Na (2.2062, 
2.2090 /xm), Sc (2.2058, 2.2071 /im), Ti (primarily 
2.2217 and 2.2240 /xm), and Mg (2.2814 /xm) lines 
as well as the CO t; = — 2 band. All of these fea- 
tures are consistent with a single K5 spectral type 
(Te// = 4300 K) with dwarf-like surface gravity 
(luminosity class IV/ V). 

2. YLW 15 has a low-mass central protostar 
(M ~ 0.5Mo) with effective temperature, surface 
gravity, and stellar luminosity similar to that of 
a T Tauri star despite its high bolometric lumi- 
nosity (Lboi — 10 Lq). We find that the stel- 
lar luminosity of YLW 15 is 3 Lq and its accre- 
tion luminosity is 7 Lq. Assuming the magnetic 
accretion model of Shu et al. (1994), we derive 
the disk mass accretion rate of YLW 15 to be 
Md = 2.3 X 1O"'^M0 yr"i and its stellar mass 
accretion rate to be = 1.6 x 10~^Mo yr~^. 
These rates are consistent with the mass accretion 
rate of the ambient cloud gas and the stellar mass 
of YLW 15 given the hfctimes of Class I protostars. 

3. YLW 15 is highly veiled in the near-IR, 
ffc > 3. We explore the origin of this veiling and 
evaluate the contributions of the stellar surface ac- 



cretion, circumstellar disk, and circumstellar en- 
velope components. We estimate that the circum- 
stellar disk and envelope contribute equally to the 
veiling each providing ^ 45% of the needed excess 
emission. We find a significant amount (r^ ~ 0.33) 
of veiling is provided by the accretion shocks at the 
stellar surface. Although this contribution is rela- 
tively small (~ 11%) compared to the total, it is 
significantly greater than that produced by simi- 
lar processes in typical CTTSs. We determine that 
emission from the circumstellar envelope is inad- 
equate to veil absorption features formed in the 
inner dust-free circumstellar disk. We conclude 
that there are no strong absorption lines formed 
in the warm inner disk of YLW 15. 

4. The projected rotation velocity of YLW 15 
is V smi = 50 km s~^, much faster than that of a 
typical T Tauri star in the p Oph cloud. However, 
this rotation rate does not appear to be correlated 
with the X-ray flarings and variability observed in 
multi-epoch X-ray observations of this source. We 
do find that it is plausible that YLW 15 is magnet- 
ically coupled to its accretion disk at a radius of 2 
R^, in the magnetic accretion model of Shu et al. 
(1994). This suggests - and we verify the plausi- 
bility - that stellar angular momentum decreases 
by over a factor of 2 in only ^ 10^ yr between the 
protostellar and T Tauri evolutionary phases. 

5. The near-IR spectrum of WL 6 does not 
show any obvious absorption features, and we esti- 
mate that the veiling of this YSO is > 4.6. This 
high veiling and its low luminosity {L\joi — 2Lq) 
dictate that its central protostar is likely to have 
low stellar luminosity (L* ~ O.6L0) and to be low 
mass (M ~ O.IM©). 
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